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ABSTRACT 

Aims. OO Serpentis is a deeply embedded pre-main sequence star in the Serpens NW star-forming region. The star went into outburst 
in 1995 and gradually faded afterwards. In many respects its eruption resembled the well-known FU Orionis-type (FUor) or EX 
Lupi-type (EXor) outbursts. Since very few such events have ever been documented at infrared wavelengths, our aim is to study the 
temporal evolution of OO Ser in the infrared. 

Methods. OO Ser was monitored with the Infrared Space Observatory in the 3.6-100 yum wavelength range, starting 4 months after 
peak brightness and covering a period of 20 months. Eight years later, in 2004-2006 we again observed OOSer at 2.2 and 12yum 
from the ground and complemented this dataset with archival Spitzer obsevations also from 2004. We analysed these data with special 
attention to source confusion and constructed light curves at 10 different wavelengths as well as spectral energy distributions. 
Results. The outburst caused brightening in the whole infrared regime. According to the infrared light curves, OOSer started a 
wavelength-independent fading after the peak brightness. Later the flux decay became slower but stayed practically wavelength- 
independent. The fading is still ongoing, and current fading rates indicate that OO Ser will not return to quiescent state before 201 1. 
The outburst timescale of OO Ser seems to be shorter than that of FUors, but longer than that of EXors. 

Conclusions. The outburst timescale and the moderate luminosity suggest that OO Ser is different from both FUors and EXors, and 
shows some similarities to the recently erupted young star VI 647 Ori. Based on its SED and bolometric temperature, OO Ser seems to 
be an early class I object, with an age of < 10 5 yr. As proposed by outburst models, the object is probably surrounded by an accretion 
disc and a dense envelope. This picture is also supported by the wavelength-independence of the fading. Due to the shorter outburst 
timescales, models developed for FUors can only work for OO Ser if the viscosity parameter in the circumstellar disc, a, is set to an 
order of magnitude higher value than usual for FUors. 
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1. Introduction 

OOSerpentis (a 2 ooo = 18 h 29 m 49 s .l, <5 20 oo = +01°16'20") is a 
deeply embedded pre-main sequence star in the Serpens NW 
star-forming region at a distance of 31 1 pc. In 1995 OO Ser un- 
derwent a large increase in K-band flux in less than 1 year, reach- 
ing its maximum brightness in 1995 October (Hodap p et alj 
Il996bl) . The object was not visible, even at peak brightness, in 
the J-band or at shorter wavelengths, therefore it is also known 
as Serpens Deeply Embedded Outburst Star (DEOS). In the H 
and K bands the object was observable but the emissi on is domi- 
nated by scattered light. In the K-band lHodappI ( I 1 9991) monitored 
the outburst until 1998 October, and found that after the peak 
OO Ser gradually faded at a rate faster than the typical fading 
rate of FU Orionis-type objects (FUors), but slower than that of 
EX Lupi-type stars (EXors). Its K-band spectrum (a steeply ris- 
ing, smooth continuum) also differed from both FUor and EXor 
spectra, which usually exhibit absorption or emission features. 
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* Based on observations with ISO, an ESA project with instruments 
funded by ESA member states (especially the PI countries France, 
Germany, the Netherlands and the United Kingdom) with participation 
of ISAS and NASA. 



Eruptions of pre-main sequence stars are rare events, thus a 
new outburst is always noteworthy. The eruption mechanism is 
thought to be caused by e nhanced accretion from the circumstel- 
lar disc onto the star (e.g. Hartm ann & Kenyonll9 96). The close 
link between the eruption phenomenon and the circumstellar ma- 
terial makes it crucial to document the outburst also at infrared 
wavelengths, where the circumstellar dust radiates. However, 
such observing programmes are constrained by the availability 
of active infrared satellite missions. The physical analysis of the 
phenomenon is limited due to the lack of preoutburst data. The 
eruption of OOSer in 1995 provided a unique opportunity to 
collect such a dataset and carry out a multiwavelength infrared 
study of the whole outburst event for the first time. 

Triggered by the news on the outburst of OO Ser we acti- 
vated a Target of O pportunity programm e on the Infrared Space 
Observatory (ISO. iKessler et all 11996). The infrared monitor- 
ing started 4 months after the maximum brightness of OO Ser, 
and continued for 20 months. The ISO-SWS measurements 
from t his programme were published separately bv lLarsson et al.l 
(2000). They found that OOSer changed its infrared fluxes in 
the 2-45 pm range and estimated an extinction of A y « 20 mag 
from the optical depth of the 10 pm silicate absorption feature. 
In an independent programme the Serpens core was surveyed by 
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the ISOCAM instrument prov iding 6.7 and 14. 3 fim photometry 
on OO Ser for a single epoch ( Kaas et al. 2004]). 

In this paper we analyse our ISOPHOT and ISOCAM ob- 
servations from the monitoring programme. In addition to ISO 
measurements, we observed OO Ser from the ground at 2.2 fim 
in 2004 and 2006 as well as at 12 fim in 2004. We complemented 
this database with archival Spitzer data also from 2004, as well 
as with previously published measurements on OO Ser from the 
literature. 



2. Observations and data reduction 

2. 1 . ISOPHOT observations 

ISOPHOT, the photometer on-board ISO (iLemke et al.lll996l) . 
carried out multi-filter photometry with 9 different filters in 
the 3.6-200 /mi wavelength range and spectrophotometry in 
the 2.47-1 1.62 /mi range, at 8 different epochs between 1996 
February and 1997 September. Table Q] shows the log of the 
observations. In most cases small raster maps were obtained 
except in 1997 September, when at 25 fim and shortwards the 
source and background positions were observed separately. The 
typical integration time was 64 s. Aperture sizes varied accord- 
ing to the filters: 1378 at 3.6 /mi, 18" at 4.8 /an, 52" at 12, 15 
and 25 /mi. At 60 and 100 /an the CI 00 camera (3x3 pixel, 
43"x43" per pixel), while at 170 and 200 /on the C200 camera 
(2x2 pixel, 89"x89" per pixel) was utilised. Some far-infrared 
observations were performed in the PHT 32 mode, which pro- 
vide d higher spatial reso lution (for a description of this mode 
see iTuffs & Gabrieill2003h . 

2.2. ISOPHOT data processing 

Standard processing and absolute flux calibration. The 

data reduction was performed using the ISOPH OT Interactiv e 
Analysis Software Package V10.0 (PI A, iGabriel et al l 1 19971) . 
The integration ramps were corrected for non-linearities. Cosmic 
particle hits were removed using the two thresholds deglitching 
method, and signal values were derived by fitting a first order 
polynomial to each ramp. The signals were transformed to a 
standard reset interval, then an orbital dependent dark current 
was subtracted and cosmic ray hits were again checked. In case 
the signal did not fully stabilise during the measurement time 
due to detector transients, only the last part of the data stream 
was used. This was found mainly in observations with the 12 
and 25 /mi filters, while at other wavelengths the measurements 
showed sufficient stability. The absolute flux calibration of mea- 
surements at 25 /mi and shortwards was done by adopting the 
default responsivity of the detector. The ISOPHOT Handbook 
dLaureiis et al.| [2003) gives a typical absolute flux uncertainty of 
40% at 3.6, 4.8, 12 and 15 /mi, and 10% at 25 /mi for ON/OFF 
staring measurements. Since our small raster maps enabled bet- 
ter background subtraction than the ON/OFF mode, the error of 
our measurements is probably smaller. Nevertheless, we adopted 
the abovementioned values as conservative estimates for the ab- 
solute flux uncertainty (see Col. 9 of Table[TJ. 

Relative flux determination. Due to the monitoring strategy, 
the observations were carried out with identical instrument 
setup at most epochs. This enabled us to determine relative 
flux variations at a certain wavelength more accurately than the 
standard processing by applying a new method developed at 
Konkoly Observatory especially for the fine relative calibration 



of ISOPHOT data (lJuhasz et al.ll2007l) . The algorithm compares 
the detector transient curves in a measurement sequence with 
the corresponding curves of a reference day (in our case 1996 
October 24) and determines a scaling factor between the sig- 
nals of the two epochs. Evaluating all epochs, one obtains a light 
curve normalized to the reference day. The final fluxes are com- 
puted as the product of the scaling factors and the absolute flux 
level of the reference day. The method provides relative uncer- 
tainties, which correspond to the formal errors of the scaling fac- 
tors. In our case, the relative uncertainties were usually around 
or below 10% (see Col. 10 of TableQ]). The absolute flux level of 
the whole light curve is determined by the absolute calibration 
of the reference flux. 



Far-infrared maps. The far-infrared observations at 60 and 
100 fim were processed with PIA in a standard way. PHT32 ob- 
servations at 170 and 200 /im were reduced using a dedicated 
software pa ckage (P32TOOLS) de veloped at MPI Kernphysik in 
Heidelberg (Tuffs & Gabriel 2003). This tool provides adequate 
correction for transients in PHT32 measurements. Absolute cali- 
bration was done by comparing the source flux with the on-board 
fine calibration source. At 60 and 100 /mi each 9x3 map was 
flat-fielded using the first raster step (3x3 pixel) as the sky flat 
position. Then the sum of two point spread functions, centred 
on the positions of OO Ser and a nearby submillimetre source, 
SMM 1 (see Sec. 12.71 ). was fitted to the brightness distribution 
on the map using the ISOPHOT measured footprint maps. A 
similar fitting procedure was applied to the PHT32 maps. The 
errors given in Table Q] represent the quadratic sum of the for- 
mal uncertainties of the fits (in the range of 3-15%) and the 
photometric calibratio n uncertainty of the d etector given in the 
ISOPHOT Handbook dLaureiis et alj|2003l) . We note that with 
this technique OO Ser could be separated from SMM 1 but not 
from other closer nearby sources, which may contaminate the 
flux of OOSer. 9x3 maps at 60 and 100 fim were obtained 
at 7 different epochs, but the first two measurements (in 1996 
February and April) were executed at late orbital phases during 
ISO's orbit. This caused a large uncertainty in these measure- 
ments, thus we decided not to present them. In the PHT32 over- 
sampled maps more than one detector pixel observed OO Ser; 
their independent photometric results were combined with a ro- 
bust averaging technique described by Abraham et alj d2004a). 



ISOPHOT-S. Spectrophotometric observations, also in the form 
of 3 x 1 small raster maps centred on OO Ser, were obtained 
with the ISOPHOT-S subinstrument. The processing of these 
data deviated from the standard scheme implemented in PIA. 
ISOPHOT-S has a successful dynamic calibration for staring ob- 
servations but in PIA this method is not applicable for rasters, 
leading to uncorrected transients and consequently reduced pho- 
tometric accuracy. Since the contrast between OO Ser and the 
background level is relatively low, we treated these rasters as 
long staring observations, and applied the dynamic calibration 
by modifying some PIA routines. The resulting photometric un- 
certainties are in the range of 5-10%. 



Colour corrections. Colour corrections were applied to each 
broadband photometric measurement by convolving the ob- 
served spectral energy distribution (SED) at a certain epoch with 
the ISOPHOT filter profile in an iterative way. The result of the 
ISOPHOT photometry is presented in TableQ] 
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Table 1. Log of ISOPHOT and ISOCAM observations. SP stands for spectrophotometry. "Map size" indicates the sizes of the 
final maps. A denotes the increment between adjacent pixel positions in the map (mapping) or the separation between source and 
background positions (ON/OFF). At certain wavelengths the beam contained nearby sources; for a detailed discussion see Sec. 12.71 
All fluxes are colour corrected. The last two columns give the uncertainties of the absolute and the relative flux calibration (Sec. 12. 2b . 



2.3. ISOCAM observations and data processing 

The 6.7 and 14.3 n m photometry was obtained with the 
ISOCAM instrument dCesarskv et al.ll 19961) on 1997 September 
22 (see Table [TJ. The data were reduc ed with the CAM 
Interactive Analysis Software V5.0 (CIA, lOtt et all [19971) . A 
dark current correction was applied following the 'VilSpa' 
method. Glitches were removed using the 'multiresolution me- 
dian transform'. This efficiently removes glitches based on the 
fact that glitches in general are much shorter than the signature of 



a real source. The data were stabilised using the Fouks-Schubert 
model. Next the individual frames were averaged to four mean 
images, one for each sky position. After this, a second deglitch- 
ing was applied based on the overlapping projected sky positions 
of these mean images. Finally, the images were flat-fielded and 
combined into the final mosaic. The pixel values were converted 
to mJy/arcsec 2 using the tabulated conversion factors available 
in CIA. 
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Photometry of the sources was obtained using the IRAF tool 
'xphot'. We used a small aperture photometry of 4'.' 5 and 6" ra- 
dius for the 6.7 and 14.3 /mi mosaic, respectively. These values 
were then corrected for the flux falling in the wings of the point 
spread function outside the chosen aperture. The small aperture 
was selected because of a nearby source (V370Ser at a distance 
of ~ 11") which otherwise would have contributed to the mea- 
sured flux. The results are shown in Table Q] The insignificant 
background level of 1 and 3 mJy/arcsec 2 does not contribute to 
the measurement error, which should be dominated by the abso- 
lute flux calibration uncertainty of 3.3% at 6.7 fim and 4.8% at 
14.3 fim (iBlommaert et al.ll2003l) . 

2.4. Near-infrared observations and data processing 

Ks -band images were obtained by the LIRIS instrument on 2004 
June 1 1 (as part of the LIRIS Guaranteed Time programme) and 
on 2006 May 6. LIRIS is an infrared came ra/spectrograph built 
at the Instituto de Astroffsica de Canarias iAcosta-Pulido et all 
120031; iManchado et alj |2004|) . and is mounted on the 4.2 m 
William Herschel Telescope at the Observatorio del Roque de los 
Muchachos (Canary I slands). The detector u sed is a 1024 x 1024 
HAWAII-1 detector (Hoda pp et alj Il996ah . which provides a 
plate scale of 07250/pixel and a total area of 4^27x4:27. In 
2004 a three point dither pattern was used with a total expo- 
sure time of 140 s, while in 2006 a five-point dither pattern was 
used and the total exposure time was 250 s. The images were re- 
duced using the IRAF package 'liris_ql' developed by the LIRIS 
team within the IRAF environment. The reduction steps consist 
of flat-fielding, sky subtraction and image co-addition. 

In Fig.Q](left) a part of our K s -band LIRIS image from 2004 
can be seen showing the surroundings of OO Ser. The eruptive 
star was clearly detected in Ks and a faint nebulosity around the 
star can be seen as well (see also Sec. 13. II ). Because of this neb- 
ulosity, photometry should be done carefull y. In order to be able 
to compare our measurements with those of Hodapp ( 1999), we 
used an aperture diameter of 1 1'.'3, the same as lHodappl(ll999h . 
Conversion from instrumental magnitudes to real magnitudes 
was done using the 2MASS J-Ky colour versus AKs (the differ- 
ence between the instrumental and the 2MASS Ky magnitudes) 
relationship for 20 comparison stars in the field. The resulting 
values are K s = 14.3+0.2 mag in 2004 and K s = 14.1 ±0.2 
mag in 2006. We note that the colours of the 20 comparison stars 
used in this calibration procedure covered a large enough range 
to include also the colour of OOSer itself. Similar p hotome- 
try wa s derived for some nearby sources (similarly to Hodapp 
d 19991) . we used a larger, 1 1'/3 diameter aperture for the two ex- 
tended objects: OOSer and V371 Ser, and a smaller 1'.'9 diam- 
eter aperture with the corresponding aperture correction for the 
point-like objects: V370 Ser, EC 38 and SMM 9); the results can 
be seen in Table [2] 

A K-band spectrum of OOSer was obtained using LIRIS 
on 2006 May 6. The observation was performed following an 
ABBA telescope nodding pattern. The total exposure time was 
2400 s, split in 4 individual exposures of 600 s. In order to re- 
duce the readout noise, the measurements were done using mul- 
tiple correlated readout mode, with 4 readouts before and after 
the integration. We used a slit width of 1" and a medium resolu- 
tion sapphire grism which yielded a spectral resolution of 2500. 
The wavelength calibration was provided by observations of an 
Argon lamp available in the calibration unit at the A&G box 
of the telescope. In order to obtain the telluric correction, the 
nearby A0V star HIP 90123 was observed with the same config- 
uration as the object. The data were reduced and calibrated using 
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Fig. 2. Normalised K-band spectra of OO Ser. Con tinuous line: 
this work, dashed and dotted lines: iHodappI d!999h . For clarity, 
we multiplied the spectrum from 1998 by 10 and the spectrum 
from 1995 by 100. 



the package iiris_ql'. Consecutive pairs of AB two-dimensional 
images were subtracted to remove the sky background, then the 
resulting images were wavelength calibrated and flat-fielded be- 
fore registering and coadding all frames to provide the final com- 
bined spectrum. A one dimensional spectrum was extracted with 
the IRAF 'apall' task. The extracted spectrum was divided by a 
composite to eliminate telluric contamination. This composite 
spectrum was generated from the observed spectrum of the cal- 
ibration star, divided by a stellar model and convolved to our 
spectral resolution. The resulting normalised K-band spectrum 
is shown in Fig. [2] 

2.5. Mid-infrared observations and data processing 

OOSer was imaged with TIMMI2 mounted on the ESO 3.6m 
telescope at La Silla on 2004 October 21, under clear and stable 
conditions. The N11.9-OCLI filter was used which has a cen- 
tral wavelength of 11.66/vm and a FWHM of 1.1 6 /mi. The to- 
tal integration time was 4 minutes. Both chopping and nodding 
amplitudes were 10", and the 072 pixel scale of the 240x320 
Raytheon detector was set. This resulted in a chop-nod corrected 
image with two negative and two positive beams captured on the 
detector. The four beams were used for independent determina- 
tion of the source flux and error. The resulting fluxes are 0.64 Jy 
for OOSer and 0.18 Jy for V370Ser (a nearby young star also 
visible in the field, see also Table |2). The flux calibrator was 
the photometric standard HD 96171 which was observed with 
the same set-up. Other photometric standard star observations of 
the same night were inspected to estimate the total absolute flux 
uncertainty of ~15%. 

2.6. Spitzer observations and data processing 

OO Ser was measured with the IRAC and MIPS instruments of 
the Spitzer Space Telescope on 2004 April 4 and 5, respectively, 
as part of the legacy programme 'c2d' (PI: Neal J. Evans II). The 
third data release of this programme containing enhanced data 
products and catalogues can be downloaded from the Spitzer 
websit^B As an example we plotted parts of the 8 fim IRAC and 
the 24 fim MIPS maps in Fig.Q] 



http://ssc.spitzer.caltech.edu/legacy/all.html 
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Fig.l. OOSer and its surroundings at different infrared wavelengths. V370Ser (EC 37), V371 Ser (EC 53), EC 38, SMM 9 and 
SMM 1 are also marked. In the middle panel, the spots below the brightest sources are instrumental artifacts ('bandwidth effects'). 
The white star in the right panel indicates the position of IRAS18272+01 14 as given in the IRAS Catalogue of Point Sources. 



The 'c2d' catalogue contains IRAC fluxes for OO Ser, but it 
is treated as a point source. The nebulosity seen in the K-band 
images, however, are still visible at 3.6 and 4.5 yum too. In order 
to be able to compare the IRAC fluxes with the ISO measure- 
ments, we decided to extract fluxes from the IRAC images for 
OO Ser using the same apertures as ISOPHOT used (diameter of 
13'.'8 at 3.6 fan and 18" at > 4.5 /mi) 

The current version of the 'c2d' catalogue does not contain 
MIPS fluxes yet. In order to obtain photometry in the MIPS 
bands, we downloaded the enhanced MIPS images. At 24 fim we 
selected 7 isolated stars to construct the point spread function. 
Since OO Ser was saturated, this profile was then fitted to the 
non-saturated wings of OO Ser. At 70 fim OO Ser and a nearby 
young star, V370 Ser could not be fully separated. Thus an aper- 
ture of 20" was utilized, which included both objects. The mea- 
sured flux was then distributed between OO Ser and V370 Ser in 
the same ratio as their peak brightnesses (measured in a 5" aper- 
ture at the position of the two sources). We also extracted 70 fim 
photometry for some other sources in the vicinity of OO Ser, 
since they possibly cause source confusion in IRAS and ISO 
far-infrared measurements. In these cases, a 20" aperture was 
utilized with a fixed 0.142 Jy background. 

Colour corrections were applied to each measurement for 
each source by convolving the observed SED with the IRAC and 
MIPS filter profiles in an iterative way. The resulting fluxes and 
estimated uncertainties can be seen in Table[2] 

2. 7. Source confusion 

Fig. Q] shows that several infrared and submillimetre sources 
are present in the vicinity of OOSer: V370Ser (also known 
as EC 37), V371 Ser (also known as EC 53 or SMM 5), EC 38, 
SMM 9 and SMM1. At shorter wavelengths (2.2, 3.6, 4.5 and 
5.8 fim), an extended nebulosity around OO Ser can also be seen. 
Thus, when comparing the fluxes of OO Ser measured with dif- 
ferent instruments one has to keep in mind that-to some extent- 
the nebulosity and some of the abovementioned sources may 
contribute to the observed flux at a particular wavelength. At 
2.2 yum and with the ISOPHOT/ISOCAM at 3.6, 4.8, 6.7 and 
14.3 fim the beam included OOSer only. At 12, 15 and 25 fim 
the ISOPHOT beams included OO Ser and V370 Ser. At 60, 100, 
170 and 200 fim the fluxes extracted for OO Ser include also con- 



tributions from V370Ser, V371 Ser, EC 38 and SMM 9. Recent 
infrared photometry for these sources can be seen in Table [2] It 
should be noted that fluxes presented in Table[T]contain the con- 
tributions of nearby sources as discussed above. In Sec. 13.41 we 
give a detailed description of how we corrected the ISOPHOT 
measurements for the effects of source confusion. 



3. Results 

3. 1 . Morphology of the nebula 

According to iHodappI (1999) before the outburst a triangle- 
shaped nebula west of OO Ser and a small elongated nebula east 
of OO Ser could be seen. During the outburst these - presum- 
ably reflection - neb ulae became much brighter. In this phase 
HodaB£etdJ 0996b) observed OO Ser in the L' and M' bands, 
as well as at 11.7 and 20.6 fim. They found that the nebulos- 
ity can be seen in the L' but not at longer wavelengths. Our re- 
cent Kj-band images reveal that the ne bulae still exis t and look 
very similar to the preoutburst image of Hodapp] (1 19991) (see his 
Fig. 1 left). The nebulosity around OO Ser is also visible in the 
Spitzer/IRAC maps from 2004 at 3.6,4.5 and 5.8 fim and some 
extended emission can be suspected even at 8 fim. 

3.2. Preoutburst fluxes 

In order to study the consequences of the outburst in the 
whole infrared regime, one has to compile first the SED of 
OOSer in the quiescent phase, i.e. estimate the preoutburst 
fluxes. At 2.2 /im there exis ts a preoutburst measurement from 
1994 (iHodapp et all Il996bl) . At longer wavelengths, only the 
IRAS measurements a re available from 1983 . Analysing high- 
resolution IRAS maps. lHurt & Barsonvld 19961) derived fluxes of 
0.63, 4.5, 24 and 111 Jy at 12, 25, 60 and 100 fim respectively. 
The authors claim that these values represent the total fluxes 
from a region encompassing three confused sources: OOSer, 
V371 Ser and SMM 9, and they give the one third of the above- 
mentioned values as upper limits for the brightness of OO Ser. 
The position of IRAS18272+01 14 from the IRAS Catalogue of 
Point Sources (marked by a white star in Fig.[T]right) is located 
halfway between the sources, which may indicate that there was 
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Table 2. Log of observations from 2004-2006. All fluxes are presented in mJy. Spitzer fluxes are from the third delivery of data from 
the 'c2d' legacy project, except the 70 fim data of all sources and the 3.6, 4.5, 5.8, 8 and 24 fim data of OO Ser, which were extracted 
by us from MIPS and IRAC images improved and published by the 'c2d' legacy team. All Spitzer fluxes are colour corrected. 



no dominant source, but all sources had comparable contribu- 
tions to the IRAS flux. 

With the help of Spitzer/MIPS measurem ents, it is possi- 
ble to check the validity of the assumption of iHurt & Barsonv 
( 1996) via estimating the preoutburst fluxes of OO Ser at 25 and 
60 fim. Assuming that EC 38 and SMM 9 have non-variable far- 
infrared fluxes, we subtracted the contribution of these sources 
from the IRAS values cited above. Another nearby source, 
V37 1 Ser, exhibits near-infrared variabilit y of » 1 .5 mag and 
shares many characteristics with EXors (Hodapp 1999). At far- 
infrared wavelengths, however, eruptive young stars do not typ- 
ically show significant flux changes (Abraha m et al.l [2004a). 
Thus we assumed that V37 1 Ser is also non-variable at 25 and 
60/im, at least within our measurement uncertainties, and we 
subtracted its flux from the IRAS values. In practice, we es- 
timated 25 fim fluxes for these nearby sources via interpolat- 
ing from the 24 fim MIPS values in Table [2] a nd subtracted the 
sum o f these (3.6 Jy) from the value given by IHurt & Barsonvl 
d 19961) (4.5 Jy). The result is 0.9 Jy, which is indeed of the or- 
der of one third of 4.5 Jy. The result at 60 fim is also roughly 
consistent with the one third value. Due to the lack of recent 
12 or 100 fim data for all nearby sources, the same test cannot 
be done at these wavelengths. Therefore, for homogeneity, at 
all four IRAS wavelengths we adopted as preoutburst fluxes the 
one third values within a factor of 2 uncertainty: 0.21^^ Jy at 
12 urn, 1 .5+^ Jy at 25 fim, 8+* Jy at 60 fim, 37+^ Jy at 100 fim. 
We note that OO Ser is the first eruptive young star where pre- 
outburst fluxes are available in the whole infrared wavelength 
regime. Similar preoutburst data can be found only for one other 
young eruptive star, V1647 0ri (Abrah am et al.l l2004b). 



3.3. Spectral energy distribution 

In Fig. [3] the SED of OOSer in the high (outburst) state from 
1996 September is plotted with filled symbols. We attempted to 
correct for the effects of source confusion as described in de- 
tails in Sec. 13.41 Thus, the resulting SED in Fig.[3]represents the 
flux from OO Ser alone at < 60 fim. Due to high extinction the 
star is invisible at optical wavelengths and is very faint in the 
near-infrared regime. In the mid-infrared (3-25 fim) the SED is 
rising towards longer wavelengths. This wavelength regime is 
zoomed in Fig.|H showing a broad silicate absorption feature at 
9.7 fim, as well as several ice features, which also indicate high 
extinction. Measuring the optical depth of the 9.7 fim feature, 
and convertin g it to visual extinction assuming Ay/ATgj = 18.5 
(|Prainel2003l) . we obtained Ay =42 + 5 mag. This va lue is twice 
as large as the one measured bv lLarsson et alj (|2000); the differ- 
ence is probably related to the low signal-to-noise of his SWS 
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Fig. 3. Spectral energy distribution of OO Ser. Squares: preout- 
burst fluxes measured with IRAS in 1983 dHurt & Barsonv 
1 19961). and K-ban d photometry from 1994 August, 
dHodapp et al.1 Il996bl) : Dots and line: outburst fluxes from 
1996 September, measured with ISOPHOT; Stars: current fluxes 
from 2004, measured with LIRIS, TIMMI2 and Spitzer. Error 
bars smaller than the symbol size are not plotted. ISOPHOT 
beams at 100 and 200 fim contained nearby sources; for a 
detailed discussion see Sec. l2.7l and l3.3l 



spectra around the 9.7 fim absorption feature. The measured high 
extinction indicates that OO Ser is more deeply embedded than 
most FUors. At far infrared wavelengths (60-200 fim) the SED 
is flat, but it should be noted that data points at 100 and 200 fim 
are contaminated by source confusion (see Sec. 12.7b . 

In Fig.|3]the preoutburst SED (see Sec. I3.2t and the present 
SED (based on measurements from 2004) are also plotted. They 
will be discussed in Sec. 13.4. 2l and l3.4.6l 

3.4. Light curves 

ISOPHOT, ISO CAM, LIRIS, TIMMI2, Spitzer, IRAS 
Hurt & Barsoiiyl 1 19961) and K-band measurements 
Hoda pp etal J Il996bt iHodappI 1 19991) were combined to 
construct the light curves of OOSer at different wavelengths 
between 1995 and 2006. 

For the subsequent light curve analyses, we attempted to cor- 
rect for the effects of source confusion. One possibility would be 
to smooth all data to the same resolution, in most cases defined 
by the ISOPHOT aperture. In doing so, however, several addi- 
tional sources would be included in the beam, falsifying the fad- 
ing rate calculations for OOSer. Instead, we decided to use the 
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Fig. 4. Ice features in the spectrum of OO Ser. Data below 12 fim 
were taken by ISOPHOT-S. Solid line indicates data from 1996 
October, dashed dot line from 1997 April. Data above 12 /mi 
were taken by t he ISO-SWS instrume nt (highly processed data 
products from Frieswi ik et alj [2004h . In order to obtain high 
signal-to-noise ratio, ISO-SWS spectra from 1997 March, April 
and September were averaged and a 0.3 fim wide moving median 
window was applied. 



higher spatial resolution images to correct for the contribution of 
additional, unrelated sources in the large ISOPHOT beams. 

At 12 /mi TIMMI2 could resolve OOSer and V370Ser. 
Assuming that V370 Ser is not variable at this wavelength, we 
subtracted its contribution of 0.18Jy (see Table |2]l from each 
ISOPHOT 12 fim measurement. The MIPS camera of Spitzer at 
24 /mi could separate OOSer from V370Ser giving a flux of 
1.56 Jy for the latter source, which we interpolated to 25 //m 
(1.71 Jy) and subtracted from the ISOPHOT points at 25 //m. 
At 60 //m, MIPS 70 fim measurements could be utilized. Using 
the SEDs presented in Tabled we interpolated 60 fim fluxes for 
V370Ser, V371 Ser, EC 38 and SMM9, and subtracted the sum 
of these values (24 Jy) from the ISOPHOT 60 //m data points. 
Due to the lack of 100 fim fluxes for the nearby sources, a similar 
correction was not possible in the case of the 100 fim ISOPHOT 
light curve. 

In Fig. [5] six representative light curves between 2.2 and 
60 fim are shown. All these light curves are corrected for source 
confusion and represent the brightness evolution of OOSer 
alone. In the following, we describe these light curves in detail. 

3.4.1. The outburst history at2.2/zm 

FUor and EXor ourbursts were historically monitored at optical 
wavelengths. Since OO Ser is invisible in the optical, Ks is the 
shortest available band where the outburst could have been fol- 
lowed. The top left panel of Fig.[5]shows the 2.2 /mi light curve 
of OO Ser. The star brightened by 4.6 mag between 1994 August 
and 1995 July. After reaching peak brightness, it started an ap- 
proximately exponential fading with a rate of 1 .00 mag/year in 
the first 350 days and 0.34 mag/y ear afterward s, as the data be- 
tween 1995 and 1999 indicate dHodappI [19991) . The change in 
fading rate divides this period of the outburst into a first and 
a second part. Our measurements from 2004-2006 prove that 
the fading continued, although at a slightly different rate, rep- 
resenting a third part of the fading. Comparing the preoutburst 
flux with the new observations, one can conclude that at 2.2 fim 
OO Ser is still above the preoutburst flux level. 
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Fig. 6. Initial rise: flux ratio of peak brightness to preoutburst 
brightness at each wavelength. (See discussion in Sec. 13.4.21 ) 



In the following, we describe the lightcurves of OO Ser at 
longer wavelengths, following the abovementioned division: the 
initial rise (until peak brightness); the first part of the fading (un- 
til mid-1996); the second part of the fading (until mid-1997); 
and the third part of the fading (until 2004). 

3.4.2. The initial rise (until peak brightness) 

Comparison of the preoutburst fluxes with the SED in outburst 
(Fig.0 shows that the eruption caused brightening in the whole 
near- to far-infrared spectrum. Though the fluxes at 100 /mi are 
contaminated by nearby sources, there is a flux change at this 
wavelength, too. The shape of the SED changed significantly: 
the outburst SED is flatter. 

At 2.2 /mi the light curve reached its peak in 1995 October. 
At 3.6, 4.8 and 12 /mi, the exact date of the peak is not known, 
but it happened not later than 1995 October, thus it was probably 
simultaneous with the K-band peak. At 25 /im, however, the peak 
took place in 1996 April, i.e. some 200 days later than at 2.2 fim. 
The shape of the peak at 25 fim is also different, it is broader 
than at 2.2 fim, and has a triangle-like shape. At 60 and 100 fim, 
the photometry is more uncertain, thus the peak dates cannot 
be determined. Nevertheless, the time-shift at 25 fim gives a hint 
that the peak brightness happened gradually later at longer wave- 
lengths. To our knowledge, such a time-shift has not been ob- 
served for other eruptive stars. 

In order to analyse the wavelength dependence of the bright- 
ening, in Fig. [6] we plotted the ratio of the peak flux to the pre- 
outburst value at 2.2, 12, 25, 60 and 100 fim. Since at 100 fim 
both the preoutburst IRAS flux and the peak ISOPHOT flux in- 
clude contributions from several nearby sources, the ratio de- 
rived from these numbers represents a lower limit for the bright- 
ening of OO Ser itself. Figure [6] suggests that the amplitude of 
the flux increase has a characteristic wavelength dependence: 
the flux ratio is lower at longer wavelengths, with an approx- 
imately linear dependence on log A. Although the peak bright- 
ness did not occur simultaneously at different wavelengths, this 
trend can also be seen if one compares in Fig. [3] the preout- 
burst SED with the SED from 1996 September (close to peak 
brightness). We note that this result is different from the case 
of another young eruptive star, V1647 0ri, where the initial rise 
was practically wavel ength-independent (Abra ham et al.l l2004b; 
Muz erolleet ai1l2005l) . 
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Fig. 5. Light curves of OOSer at 
different wavelengths. Open circles: 
measur em ents o f iHodapp et al.l 
dl996bh : iHodaprj dl999h : Stars: 
LIRIS; Filled circles: ISOPHOT; 
Triangle: TIMMI2; Squares: Spitzer; 
Dashed lines: preoutburst fluxes 
from 1994 at 2.2yum and from 1983 
at 12, 25 and 60 /mi (for details 
see Sec. 13. 2\ . Tick marks indicate 
January 1 of the corresponding 
year. Error bars in the upper right 
corner of each panel show the typical 
uncertainty in the absolute flux level. 
Due to the homogeneous measure- 
ment strategy, the uncertainty of the 
ISOPHOT data points relative to 
each other at a certain wavelength 
are usually better than the absolute 
errors (see Sec. 12.2b . All fluxes are 
corrected for source confusion (for a 
detailed description see Sec. 12.71 and 
3.4t . thus these light curves show the 
brightness evolution of OO Ser alone. 
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Fig. 7. Fading rates of OOSer at different wavelengths, (a) 
between 1996 April and 1996 September, (b) between 1996 
October and 1997 March, (c) between 1997 April and 2004. 



3.4.3. The first part of the fading (until mid-1 996) 

During the first year after peak brightness, OO Ser exhibited sim- 
ilar exponential fading at all wavelengths. This can be quanti- 
fied by fitting a linear relationship to the first part of the light 
curves plotted logarithmically in Fig. |5J using data taken only 
after peak brightness at a given wavelength. The derived fading 



rates are displayed in Fig. 0(a). This graph shows that the fad- 
ing rate during this period was approximately 1.0 + 0.3 mag/year 
irrespectively of wavelength in the whole mid-infrared regime. 



3.4.4. The second part of the fading (until mid-1 997) 

The second part of the light curves (between 1996 October 
and 1997 September) can be characterised by fading rates of 
0.35+0.10 mag/year, similar at all wavelengths from 2.2 to 
60 jum (Fig.|7](b)). This is approximately three times slower than 
the fading during the first part. In this period ISOPHOT-S spec- 
tra were also obtained, which give the opportunity to analyse 
the fading with high wavelength resolution. In order to increase 
the signal-to-noise ratio, the fitted fading rates of ISOPHOT-S 
were binned (in the case of the short wavelength channels) or 
smoothed (in the case of the long-wavelength channels). The re- 
sulting values are overplotted with small open or filled dots in 
Fig. |7](b), respectively. These values also support a wavelength- 
independent fading, characteristic of the second part of the out- 
burst (except where the emission of OO Ser is strongly reduced 
by the silicate absorption around 10/im, introducing a large un- 
certainty in the fit). 



3.4.5. The third part of the fading (until 2004) 

New measurements from 2004 revealed that the fading moni- 
tored by ISOPHOT continued at all wavelengths after 1997. We 
calculated fading rates for the period 1997-2004. The resulting 
numbers, which are plotted in Fig. [7] (c), are even lower than 
those for the second part and they also show no wavelength de- 
pendence. The values are in the range of 0.19 ± 0.08 mag/year. 
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3.4.6. The present status of OOSer 

In Fig.[5]the preoutburst fluxes (marked by dashed lines) can be 
compared with the latest measurements from 2004-2006. This 
comparison reveals that in 2004 OO Ser was still above the pre- 
outburst state in the whole infrared wavelength regime, indicat- 
ing that the outburst had not finished yet. The same conclusion 
can be drawn from Fig. [3] if one compares the SED from 1996 
with that from 2004. 

In Fig.[5]the last K-band photometric point from 2006 seems 
to deviate from a linear extrapolation of the preceding fading, 
possibly indicating a fourth phase of the outburst with very slow 
flux change. Extrapolating the lightcurves using the rates calcu- 
lated for the third part of the fading (1997-2004), the expected 
end date of the outburst is w 201 1 from both the 2.2 /jm and the 
12 fim lightcurve. The slowing down of the fading process, how- 
ever, may delay this event even well beyond 201 1. 

3.5. K-band spectral evolution 

Figure|2]displays our K-band spectrum from 2006, together with 
the outburst spectra from 1995 and 1998 (iHodappll 1999b . All 
three spectra are taken with comparable slit widths and were 
normalised to their value at 2.4 /mi. In general, all spectra show a 
steep continuum, rising towards longer wavelengths. No individ- 
ual absorption or emission lines can be seen, even in quiescence. 
The overall shape of the three spectra is very similar, though 
there is a tendency that later spectra are less steep. 

4. Discussion 

4.1. The outburst 

OO Ser offers a unique possibility to investigate the long-term 
behaviour of an eruptive young stellar obj ect. The rise time o f 
OOSer was 8.5 months at the longest dHodapp et al.| [l996b). 
similar to the typical timescale of "fast" FUor outbursts (e.g. the 
outburst of FU Ori or VI 057 Cyg, w hich was modelled by a trig- 
gered eruption, see iBell et alll995l) . 

As discussed in Sec. l3.4.fl the fading rate abruptly decreased 
350 days after the maximum brig htness. A similar event was 
observed in the case of V 1 057 Cyg dKenvon & Hartmandfl99"I 
too, but there the transition happened later, a bout 1300 days after 
the m aximum. Using the light curves of iKenvon & Hartmannl 
( 1991) we computed fading rates for VI 057 Cyg and compared 
them with corresponding values of OO Ser. We found that both 
before and after the transition, the fading of OO Ser ( 1 .00 mag/yr 
before, 0.35 mag/yr after) was significantly faster than that of 
V1057Cyg (0.09.. .0.40 mag/yr before, 0.04.. .0.14 mag/yr after, 
between 0.44 and 21 ^im). 

From the 2.2 //m light curve in Fig. [5] we can conclude that 
the object will return to the quiescent phase some time after 
2011. This implies that the duration of the outburst of OOSer 
will be at least 16 years. This timescale differs both from that of 
FUors (being several decades or a century) and that of EXors 
(being some weeks or months), suggesting that OOSer is a 
young eruptive object that differs from both FUors and EXors. 

The magnitude of the luminosity change during outburst 
is another argument in favour of OO Ser being different from 
FUors or EXors. From Fig. [3] we could estimate a bolomet- 
ric preoutburst luminosity of L pre = 4.5 + 1 .5L Q for OO Ser and 
from the ISOPHOT data we also computed luminosities for 
OO Ser for the different epochs during the outburst. Since the 
far-infrared fluxes are contaminated by nearby sources, we cal- 



culated an upper and a lower limit for the luminosity, by includ- 
ing or neglecting the 100 - 200 //m data points, respectively. 
Adopting the 1996 February values of L pea k =26 . . . 36L as a 
representative outburst luminosity range, OO Ser changed its lu- 
minosity by a factor of about 6 ... 8. Thus one may conclude 
that both the peak brightness and the amplitude of the lumi- 
nosity increase was significantly lower than the corresponding 
values of classical FUors (~ 100 L Q and a factor of 100, see 
lHartmann & Kenvonll 19961) . 

We note that there exists another star, VI 647 Ori, which 
seems to share some characteristics of OOSer. Based on 
the amplitude of its brightening, VI 647 Ori was classified 
as an intermediate-typ e object between FUors and EXors by 
iMuzerolle et alj d2005l) . Preoutbu rst and outburst lumino sities 
of V1647Qr i are L pre = 5.6L dAbraham et alj 12004b) and 
Lpeak = 44 L Q dMuzerolle et al.ll2005l) . respectively, thus its lumi- 
nosity changed by a factor of about 8, similarly to OOSer. The 
timescales of their outbursts are somewhat di fferent, because the 
eruption of V 1647 Ori was only 2 years long (Kospal et al. 2005) 
and it also produced an outburst in the 1960s dAspin et alj2 006). 

We speculate, following IHodapp et ail d!996bl) . that OO Ser 
(and in some respects VI 647 Ori) may be the representative of a 
new class of young eruptive stars ("Dee ply Embedded Outburst 
Star" or DEOS in IHodapp et al.lll996bT) . Members of this class 
may be defined by their relatively short timescales compared to 
FUors, possibly recurrent outbursts, modest increase in bolomet- 
ric luminosity and accretion rate, and an evolutionary state ear- 
lier than that of typical FUors or EXors (see Sec. 14.21 ). 

4.2. The evolutionary stage of OO Ser 

Base d on optical to milli metre measurements available at that 
time, IHodapp et all d!996bl) claimed that the SED of OO £>er ap- 
peared consistent with that of a class I source and assumed its 
age to be 10 5 yr, indicating that OO Ser is in a very early evo- 
lutionary phase. The new data presented in this paper make it 
possible to reestablish th e evolutionary stag e of the source. We 
followed the method of IChen et alj 0995) and computed the 
bolometric temperature 7boi according to their Eq. (1) for the 
preoutburst, outburst and present SEDs. Due to the uncertainty 
of the far-infrared data points (see Sec. l2.7l about source confu- 
sion), we calculated 7boi in two different ways, with and with- 
out data points affected by source confusion (60 < A < 800 //m). 
The resulting values are in the range of 50 - 120 K for all three 
SEDs. Within this interval, the bolometric temperature slightly 
increased, and the luminosity changed by a factor of 7 during the 
outburst. We compared these values with the distribution of cor- 
responding values among young st ellar objects in the Taurus and 
p Ophiuchus star forming regions dChen et aTl [l995). From this 
check we can conclude that OO Ser seems to be an early class I 
object, and its age is < 10 5 yr. 

4.3. Structure of the circumstellar environment 

The circumstellar environment of young eruptive stars is usually 
modelled with a flat or flared accretion disc and an extended in- 
fallin g envelope (e.g. IKenvon & Hartm ann 199lt iTurner et al.l 
119971) . In these models the emission of the central source (the 
star and the innermost part of the accretion disc) dominates the 
emission at optical and near-infrared wavelengths. Between 3 
and 10 /mi the origin of the emission is the release of accretion 
energy in the disc, and also starlight reprocessed in the same part 
of the disc and also in the envelope. The emission at A > 10/im is 
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starlight reprocessed in the envelope. The outburst occurs when, 
due to thermal instability in the inner edge of the disc, the ac- 
cretion rate suddenly increases. After peak brightness, the ac- 
cretion slowly relaxes to its quiescent value. The decreasing ac- 
cretion rate causes the fading of the central source and conse- 
quently leads to the simultaneous fading of the reprocessing en- 
velope. Thus, this model predicts a wavelength-independent fad- 
ing of the source in the whole optical to mid-infrared regime 
(see e.g. V1057Cvg inlXbraham et al.l l2004al and V1647 0ri in 
lAcosta-Pulido et alJl2007F ~ 

The circumstellar environment of OO Ser probably shares 
many properties with the abovementioned models, except that 
the unusually high extinction indicates a larger and/or denser en- 
velope. As it can be seen in Fig. [7] the fading of OO Ser was in- 
deed wavelength-independent in the whole near- to mid-infrared 
wavelength regime, in agreement with the model predictions. 
This, together with the overall shape of its SED and the ice fea- 
tures in its mid-infrared spectrum indicates that the circumstellar 
structure of OO Ser is similar to those of other young eruptive 
stars, i.e. possesses a circumstellar accretion disc and is embed- 
ded in a dense circumstellar envelope. 

4.4. Viscosity 



iBell & Lid (11994 modelled FUor outbursts as self-regulated 
accretion events in protostellar accretion discs. In this model 
the risetime of the outburst, the subsequent high state and the 
time between successive outbursts are dependent on a, t he vis - 
cosity parameter in the model of IShakura & Sunvaevl (11 973). 
Fitting their model to the observed FUor timescales, [Bell & Lin 
( 1 19941) derived a c = 1CT 4 for the a value in the cool, neutral state, 
and a/, = 10~ 3 in the hot, ionised state. Sinc e in many r e spects 
OO Ser is similar to FUors, the model of IBell & Linl d 19941) 
might be applicable, although the timescales of the OO Ser out- 
burst are remarkably shorter than that of FUors. Thus, applying 
this model to OO Ser r equire s different parameters than those 
for FUors. IBell & Linl d 19941) give in their Table 2 the depen- 
dence of different timescales on a, from which we can estimate 
ah = 10~ 2 for OO Ser. This is one order of magnitude higher than 
the usual value for FUors, which implies that OO Ser may dif- 
fer from classical FUors in a way that its disc has different, one 
order of magnitude higher viscosity. 

5. Summary and conclusions 

In this paper we presented an infrared monitoring programme on 
OO Ser, a deeply embedded young eruptive star in the Serpens 
NW star-forming region. OO Ser went into outburst in 1995 and 
has been gradually fading since then. Our infrared photometric 
data obtained between 1996 and 2006 revealed that the fading 
of the source is still ongoing in the whole infrared wavelength 
regime, and that OO Ser will probably not return to quiescent 
state before 2011. The flux decay has become slower since the 
outburst peak and has been practically wavelength-independent. 
From these results we draw the following conclusions: 

- The outburst timescale and the moderate luminosity suggest 
that OO Ser is different from both FUors and EXors, and 
shows some similarities to the recently erupted young star 
V1647 0ri. 

- Based on its SED and bolometric temperature, OO Ser seems 
to be an early class I object, with an age of < 10 5 yr. As pro- 
posed by outburst models, the object is probably surrounded 
by an accretion disc and a dense envelope. This picture is 



also supported by the wavelength-independence of the fad- 
ing. 

- Due to the shorter timescales, outburst models developed for 
FUors can only work for OO Ser if the viscosity parameter 
in the circumstellar disc, a, is set to an order of magnitude 
higher value than usual for FUors. 
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